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We analyse the recent experiments investigating the low-energy physics of three lithium 6 
atoms in three different internal states with resonant two-body scattering lengths. All observed 
features are qualitatively consistent with the expected Efimov effect, i.e. the effective universal 
three-body attraction that arises for large values of the scattering lengths. However, we find that a 
quantitative description at negative energy requires non-universal two- and three-body corrections 
due to presently unknown behaviour at short distance. An attempt to implement these corrections 
is made through energy-dependent parameters that are fitted to the experimental data. 

Nous analysons les recentes experiences concernant la physique a basse energie de trois atomes 
de lithium 6 dans des etats internes differents avec des longueurs de diffusion resonantes. Toutes 
les observations s'expliquent qualitativement par l'effet d'Efimov, c'est-a-dire V attraction universelle 
effective a trois corps qui survient aux grandes valeurs des longueurs de diffusion. En revanche, 
il apparait qu'une comprehension quantitative necessite des corrections non-universelles a deux et 
trois corps provenant du comportement inconnu a courte distance. Nous proposons d'implementer 
ces corrections par le biais de parametres dependant de V energie et ajustes a V experience. 

In 1970, Vitaly Efimov predicted a peculiar effect [!Q, a universal effective attraction occurring between three 
quantum particles with short-range interactions, whenever the two-body interactions are resonant, that is to say 
their respective two-body scattering lengths are much larger than the typical short range of those interactions. 
. Interestingly, this attraction can bind the three particles to form trimers irrespective of the sign of the 
scattering lengths, i.e. whether the two-body interaction supports a two-body bound state of similar energy or 
not. In particular, there is a domain of low energies and large scattering lengths where the three-body physics is 
dominated by the Efimov attraction. Because this attraction scales as the inverse square of the mean distance 
between the three particles, as does the kinetic energy, the physics in that region is invariant under discrete 
scale transformations. As a result, for an infinite scattering length, an infinite set of trimers accumulate at 
zero-energy, and all these trimers in that region have the same structure described by only a few parameters. 
For this reason, it is referred to as the (three-body) universal region [3J. 

Within the last year, the Efimov effect has received a wealth of experimental confirmations from the ultracold 
atom community, for several kinds of atomic species [3H9]. Evidence of Efimov trimer states appeared as 
strong enhancement of inelastic collisions due to the presence of an Efimov trimer just below the collisional 
threshold. These successes were made possible thanks to the existence of magnetic-field induced Fano-Feshbach 
resonances [10] . which enable experimentalists to adjust the two-body scattering length to very large values. 
The case of lithium 6 atoms O El [3 [12j [13] is particularly striking for several reasons. First of all, unlike the 
species used in other experiments, lithium 6 atoms are fermions. Since the Efimov attraction does not occur 
for identical fermions because of the Pauli exclusion principle, they have to be prepared in three distinguishable 
states 1, 2 and 3. As a result, unlike identical bosons, this system is described by three different scattering 
lengths, and features three dimers with different binding energies. Fortunately there exists a Feshbach resonance 
for each pair in the same range of magnetic field, enabling one to have large scattering lengths for all pairs at 
the same time and the Efimov effect to manifest itself in various ways |llj . In that region, the atoms can pair 
to form three different kinds of dimers (12), (23), and (31). This creates the possibility of chemical exchange 
reactions in the universal region, such as (12)+3— >-(13)+2 where the energies intersect [Mj. Finally, because 
those dimers are comparatively stable, it is possible to study the connection between dimers and trimers at 
negative energy. The three-component lithium 6 gas is the first system where a dimer and an atom could be 
associated into one of the trimers, enabling direct trimer spectroscopy |15j . 

In this paper, we examine in detail how to theoretically reproduce the experimental observations of the Efimov 
effect in lithium 6. First, we describe the two-body physics, explaining how to model the Feshbach resonance 
occuring for each pair of atoms. Then, we explain how to deal with the three-body problem, making use of 
some approximations. We finally apply these models to analyse the experimental observations. 

I. TWO-BODY PHYSICS: THE FESHBACH RESONANCES 

A lithium 6 atom has a nuclear spin i — 1 and a valence electron spin s = 1/2, and their coupling leads to 6 
different hyperfine states in the electronic ground state. These states can be separated by applying a magnetic 
field, using the Zeeman effect. They are labelled from 1 to 6 in order of increasing energy -E^i,...^- Since 
lithium 6 atoms are fermions, we consider all the possible antisymmetrised pairs {i, j} of these states as a basis 
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of diatomic channels. The interaction between two atoms A and B depends on the arrangement of the total 
electronic spin S = sa + sb into singlet (S = 0) or triplet (S — 1) states. It can be written as: 

V(r) = Vo{r)P + Vi(r)P 1 , (1) 

where r is the distance between the two atoms, P<?=o,i are the projectors onto the singlet and triplet states, 
and Vs = o,i( r ) are the respective singlet and triplet potentials. Because the interaction is non-diagonal in the 
diatomic channel basis, it leads to a set of coupled Schrodinger equations describing the relative motion of a pair 
of atoms. As a result, a bare scattering state in one channel can be coupled to a bare bound state in another 
channel. Since these states have different magnetic moments, it is possible to change their energy difference by 
applying an external magnetic field. As these states get closer in energy, the scattering length of the dressed 
states (solution of the full equations including the coupling) becomes larger. At some point, it diverges and 
changes sign. That phenomenon is a so-called Fano-Feshbach resonance [TU]. Thanks to the existence of such 
resonances, it is possible to make the scattering lengths very large by controlling the magnetic field. 

To quantify the effects of the resonances used in the lithum 6 experiments, we solved the coupled-channel 
equations using singlet and triplet potentials for lithium. The potentials were obtained from the combination 
of ab-initio calculations and adjustments to previous experimental data from Ref. and are currently the 
most accurate potentials for lithium 6. They were kindly provided to us by Paul S. Julienne and Eite Tiesinga. 
From the solutions, we extract important quantities characterising the interaction. For scattering states of an 
atom pair {i, j} at some positive energy E — —p 2 (m is the mass of lithium 6) just above the channel threshold 
Ei + Ej, the component in the channel {i,j} (the so-called open channel) has the following asymptotic form 

, p. _. sin or , . cos pr 

^K^r-aiM—^ (2) 

which defines the energy-dependent scattering length ay(p). This quantity is directly related to the more 
familiar scattering amplitude fij(p) = —l/(ip + l/a,ij(p)). Similarly, for states of negative energy E = ~^k 2 
below the channel threshold (physical bound states as well as unphysical states), the asymptotic form becomes 

F . sinhftr cosher 

Vw^ — — — - — (3) 

which extends the definition of <2jj (p) to negative energy. The only physical states are the bound states which 
decay exponentially at large distance. Therefore they must satisfy the condition 

aij(in) — 1/k, 

whose solutions Ky correspond to the discrete spectrum of binding energies E.^ 
At low energy, the energy-dependent scattering length can be expanded as 
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p 2 + ... (4) 



from which we can extract the zero-energy scattering length and the effective range r e ,y. Note that bound 
states with small binding energies must satisfy the universal two-body property K%j ~ 1 /aij. 

The various scattering lengths, effective ranges and binding energies of the last diatomic bound states are 
represented in Figs. [I] and [2] as a function of magnetic field for the different atom pairs {1,2}, {2,3}, and {1,3} 
relevant to the experiments. According to the universal two-body property K%j ~ 1/fly , a bound state appears 
at each resonance point where the scattering length becomes infinite, and its binding energy increases as 
h 2 /ma 2 j on the side of positive values of the scattering length. This universal binding energy, only valid for 
large positive scattering lengths, has been represented as dashed lines in Fig. [T] By comparing the scattering 
lengths to the typical range of the potentials, the van der Waals length a [TU], and the binding energies to the 
van der Waals energy h 2 /ma 2 , we can define two different universal regions A and B, which are indicated in 
Fig. [T] For lithium 6, a ~ 1.5815 nm. 

For convenience, in the rest of the paper we will note Ojfc, Kju, Ejk, etc, as a,, ftj, Ei where {i, j, k} = {1, 2, 3}. 



II. THREE-BODY PHYSICS: APPROXIMATION SCHEMES 



It would be rather difficult to solve the problem of three atoms using the two-body coupled potentials described 
in the previous section. Moreover, on top of these two-body interactions, there is also a three-body interaction 
which is presently unknown. However, since we are considering the low-energy properties of the system, a 
precise knowledge of those interactions and an exact solution of the three-body problem is not necessary. At 
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Figure 1: Top panel: scattering lengths for each pair of atoms as a function of magnetic field. Bottom panel: two-body 
and three-body energy spectra, as a function of magnetic field. The solid coloured curves correspond to the dimers' 
energies, and the dashed coloured curves are their universal limit —h 2 /rna 2 j. The dashed grey curves correspond to 
the trimers' energies based on the universal theory (with |A| = 1.165 a -1 ). The orange curves correspond to the single- 
channel contact model with both two-body and three-body non-universal corrections implemented through an energy 
dependence of a%j and A. The white areas are the universal regions where all scattering lengths are larger than 3a and 
all energies are smaller than h 2 /m(3a) 2 . 
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Figure 2: Effective range of the interaction for each pair of atoms, as a function of magnetic field. Solid curves are 
calculated from our multichannel two-body model. Dashed curves are obtained from two-channel model in Eq. 1 16 1. 



low energy, the wave function has a large de Broglie wavelength and is affected by the small-scale details of the 
interactions only through some phase shifts which can be set by choosing appropriate boundary conditions at 
short distance ! 1| . Alternatively, one may replace the interactions by pseudo-potentials with the same low-energy 
effects (causing the same phase shifts) , or construct a low-energy effective field theory [T7] with a few parameters 
(corresponding again to these phase shifts). All these approaches are essentially equivalent, and are intended to 
simplify the problem by reducing the number of coordinates, and invoking just a few parameters instead of the 
complicated knowledge of the two- and three-body interactions. Here we will use the pseudopotential method. 
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We will consider two different pseudopotcntials: a single-channel contact pseudopotential, and a two-channel 
separable pseudopotential of finite range. 



A. Single-channel contact pseudopotential 



Let us start by considering a single-channel pseudopotential as a substitute for the real interaction. The 
asymptotic behaviour of the two-body scattering state in the open channel in Eqs. ^ and ^ can be reproduced 
by considering a noninteracting wave ip with energy and imposing the Bethe-Peierls boundary condition 
ip oc 1 — or equivalently (rip(r))' /(rip(r)) > —l/a(p) which fixes the logarithmic derivative of ip at 

short distance r. This condition is equivalent to a Fermi-Huang- Yang contact pseudopotential V defined by: 



(r\VU) = 4%h a ^ 5 3 (r) ^-(rtp(r)) with p 2 = lim ■ 
m or r-fQ 



tp(r) 



(5) 



By constructing such a pseudopotential Vi(Fi) for each pair of atoms {jk}, the 3-body Schrodinger equation 
reads 









[-- 

m 







E + V 1 (r 1 ) + V 2 {f2) + Vs{r 3 ) 



V(R,r)=0, 



(6) 



where we have chosen a particular set of Jacobi coordinates (R, r) = {R\,f\) to represent the relative config- 
uration of three atoms in states 1, 2, and 3 - see Fig. [3j The pseudopotentials ensure that in all three Jacobi 
coordinate systems, the three-body wavefunction ^ should locally have the following form [15] : 



*{Ri,n) 



OiiPi) 



Xi(Ri) + 0( ri ), 



(7) 



2 -V»* mE , 3^%Xi(Ri 

with = lim ' — 



which is consistent with the Bethe-Peierls boundary condition for each pair of atoms. This introduces the three 
quantities Xi(Pi)- From Eq. (JsJ , we obtain 

jV|-V?-^)*(^f)=4 7 r £ Xi(^)^ 3 (r1). (8) 

' i=l,2,3 

The solution of this equation in Fourier space is 



where the term ^Pg is a solution of the homogeneous (free) 3-body equation. It physically corresponds to an 
incident wave of 3 free atoms (2n) 6 S 3 (P — Po)S 3 (p — po) in the case of positive energy E = ^r(jP§ + Po), or 
it must be taken to be zero in the case of states with negative energy E, i.e. states \f r with at least two bound 
atoms. From Eq. (|9| and the definition of Xi, it follows that the Xi must satisfy a set of 3 coupled equations 
known as the Skorniakov - Ter-Martirosian equations |19] , 



1 \ 1 f A , q P 2 + q 2 + Pq- mE/ti 2 . , . , ^ f d 3 p ~ , -. _ 

tp j MP) - - 1 m p2 - \ 2 _ p q q _ J /h2 ixM + *.<*)) - jj^MP,P) do) 



with the relative momentum 7p = «/ 1 P 2 — The derivation is given in Appendix A. If we take the upper 



bound A of the integral in the left-hand side of Eq. ( 10 ) to be infinite, the equations admit several solutions, 



as was noted by G. V. Skorniakov and discussed by G. S. Danilov [20]. This is because the replacement of 
real interactions by two-body contact interactions is in general not sufficient to have a well-defined 3-body 
problem. In special situations of two-body interactions with a large negative effective range, a unique solution 
is determined by the equations [T5] . Otherwise, an extra boundary condition on the wave function at very short 



Figure 3: Jacobi coordinate systems (Ri,fi). In Fourier space, we have conjugate wave vectors (Pi,Pi). 



distance where all three atoms are close to one another is necessary. Cutting off the integral at some momentum 
A |22j precisely provides such a condition, as shown in Appendix B. 

In his original papers, Vitaly Efimov considered an energy-independent scattering length a(p) = a. As we 
mentioned before, at the two-body level this approximation is valid only in the universal region of large scattering 
lengths (compared to the range of the interactions). Hence, Eqs. (10) with an energy- independent scattering 
length are equivalent to the universal theory originally proposed by Vitaly Efimov to investigate the universal 
properties of 3-body systems with large scattering lengths. Note that the energy-dependence of a(k) brings 
finite-range corrections to the universal theory, including the effective range correction [23] as can be seen from 
Eq. @. 



B. Coupled-channel separable pseudopotential 

Although the previous approach does take into account finite-range corrections, it may still look oversimplified 
because it neglects the coupled-channel nature of the real atomic interactions near Feschbach resonances. To 
account for this, one may substitute the real interaction by an effective two-channel interaction U that reproduce 
most features of the original Feshbach resonances 24, 27]. It can be written as 



U = U \o)(o\ + U c \c)(c\ + U oc \c)(o\ + U co \o)(c\ (11) 

where \o) and |c) correspond to open and closed channels, respectively. 

Let us first look at two atoms interacting with this effective potential. The two-body Hamiltonian is K + U, 
where K is the relative kinetic energy. To model a resonance with a single bound state <f> of energy E\, in the 
closed channel, we take the closed-channel Hamiltonian K + U c to be Eb\<j)) (4>\ and the channel-coupling term 
as U oc — Ul Q — W\(j)){<j)\. A two-body state \ip) can be written as a superposition of open and closed-channel 
components \ip) — J rgA A(jp) \pj \o) + B\(p) |c) , where \p) is the plane wave for two atoms with relative momentum 
p, and the Schrodinger equation (K + U — E)\tjj) =0 leads to: 

(E b — E)B +2 f<£?A(q)W<f>*(<?) =0, (12) 



(2tt) 3 

( h ^^E)A(p) +J J J A(q)U o (q,p) + BWcl ) (p) = 0. (13) 



By assuming that the open-channel potential is separable [2"4Tf2"T] . U (q,p) = ^^-Xip* (q)(j)(p), we can group 
the last two terms in the second equation, which greatly simplifies the problem. This arbitrary choice is 
permitted since the small-scale details of the interaction are unimportant as long as they correctly reproduce 
the low-energy physics, as argued before. For the same reason, we can arbitrarily choose the functional form of 
4>(p). Following [SUET], we choose for convenience a Gaussian form <j)(p) = e~2 l - bp > with range b. 

Eliminating B in Eq. (13) and solving at positive energy E = — ^a, one obtains: 



A(p) = A (p)- r J iP °l P) , (14) 
— E — is 

m 

where Aq(p) = (27r) 3 <5 3 (p — po) corresponds to an incident plane wave of momentum po, and we introduced the 
two-body T— matrix element 
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Solving for T(po,p) self-consistently using the last two formulae, we deduce that the scattering length a(p) 



takes the form 



a(p) 



E b -E 



/7F6 



pErR(pb) 



(16) 



where a 



2nh 



2 |A| 2 , and Erfi is a real function related to the standard error function Erf by Erfi(z) = Erf(iz) ji. 



By adjusting the parameters A, a, E b and b, we can construct for each pair of atoms {j, k} a pseudopotential 
Ui which reproduces the energy-dependent scattering length a,i(p) of the real interaction. 

We now consider three atoms 1, 2, 3 interacting through these pseudopotentials U\, U% and U3. The 3-body 
wave function can be written as: 



l*> = 



d 3 P d 3 p 
(2tt) 3 (2tt) 3 



A(P,p)|P) 1 |p) 23 



E 

(i,j',fe)=(l,2,3) 



d 3 P 

(2tt) 3 



Bi(P)\P) 



(17) 



where |P)j is the plane wave with momentum P for the relative motion between atom i and the centre of mass 
of the pair of atoms j and k, \p)jk is the plane wave of momentum p for the relative motion between atoms j 
and k, and \4>)jk is the closed-channel bound state for the pair of atoms j and k. Defining the quantity 



i(Q) = ~ \\ 



E, 



Ai 

3 fi 2 Qi 



4 m 



- E + ie, 



d 3 p 



3 4>*{Pi)A(Qi,pi), 



we obtain an equation similar to Eq. ^ for the open-channel component A 
A(P,p) = A Q (P,p) + 3 — 



|p2 +p 2 _ m E/h 2 + ie 



i c ^ ^ 



(18) 



(19) 



2.3 



where A is an incident plane wave (2tt) 6 5 3 (P — P )S 3 (p — p Q ) for positive energy E = ~{\Pq +Po), or zero 
for negative energy. The fa then satisfy the generalised Skorniakov - Ter-Martirosian coupled equations, 



a>i{ilp) 



+ lp)fa{P)-- 

TT 



q dq 



du— 
■ 1 



P 2 + q 2 - Pq-mE/h 2 



Pqu) 

0j(Q) 



fa(q)) 



d 3 p 
(2?r)3 



(10) 



Note that these equations are very similar to Eqs. 
in the energy-dependent scattering length aj, except for the terms 



In particular, all the two-body physics is contained 
When (f> — > 1, i.e. when the range b of 
the interaction goes to zero, we formally retrieve the single-channel contact interaction equations (10). This 
indicates that apart from the terms 4> there is little difference between the two approaches. Here, the presence 
of cj) with nonzero range b effectively cuts off the integral at high momenta, so that we can safely take A — > 00. 
In fact, the range b of the pseudo-potential plays the role of A in Eqs. (10), i.e. it characterises the 3-body 
behaviour at short distance. 

It might seem unreasonable to choose 6, which is determined by the two-body interaction only, in order to 
characterise a three-body property. In general, the short-range three-body behaviour should also depend on 
a three-body interaction between the atoms. To be consistent, a three-body interaction should be added in 
the problem. Another approach, taken in Ref. [27], is to regard b not as a parameter describing the real two- 
body interaction (for example its effective range), but as a free two-body parameter that is adjusted to set the 
combined effects of the real two-body and three-body interactions. 



C. Models 



To analyse the experimental results, we will consider three different three-body models: 



The universal model is given by Eqs. ( 10 ) with energy-independent scattering lengths et^ that are obtained 
from the two-body calculation of section I, and A is a free parameter. 
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The single- channel contact model is given by Eqs. (10 1 with energy-dependent scattering lengths di(p) 
given by Eqs. (|16[), where we set the parameters A; a*, and 6 so as to reproduce the two-body 

quantities calculated in Section I, namely dj, r e ^ and Kj, as a function of magnetic field. The reason why 
we use this two-body model instead of the real Oj (p) is that it is numerically intractable to calculate the real 
di(p) at negative energy due to the divergence of sinh and cosh at large distance in Eq. Nevertheless, 
this simple two-channel model should be very close to the real Oj (p) in the low-energy domain that we are 
interested in. 



The two-channel separable model is given by Eqs. (20) with the scattering lengths aj(p) given again by 
Eqs. jl6l). 



III. ANALYSIS OF THE EXPERIMENTS 



There has been two sets of experiments. The first one [3l[T2l[T3] consisted in measuring the three-body recom- 
bination rate in a gas of atoms equally distributed in the lowest three spin states. Three-body recombination is 
the process where three colliding atoms recombine to form a combination of states with lower internal energy, 
such as a deep dimer and a free atom. At the Feshbach resonance locations, the three-body recombination is 
strongly enhanced by the very large scattering length. Away from these points, other peaks were found and 
attributed to recombination enhancement by the presence of Efimov trimers at zero energy. Indeed, whenever a 
trimer state exists just below the collisional threshold of three atoms, the three atoms resonate, which increases 
their probability to be close together and recombine. In the second set of experiments [3 [7], a gas of dimers and 
atoms was prepared, and inelastic collisions by relaxation to deeper dimer states were observed. The relaxation 
rate was also found to be enhanced at two magnetic field values due to the presence of two Efimov trimers just 
below the collision threshold. Recently, one of the two trimers' energy was directly observed by association 
spectroscopy |15l 128) . These experiments thus provided some partial information about the spectrum of Efimov 
trimers at both zero and and negative energies. The general spectrum based on these results is given in Fig. [T] 



A. Experiments at zero energy (3-body recombination) 



When three atoms in different states collide and recombine, the density n.j of atoms in each state i decreases 
according to the rate equation: 



= -KjecniUjUk, 



(21) 



where -Krec is the recombination coefficient. By measuring the variation of the number of atoms, and taking 
into account other kinds of loss, it is possible to extract the recombination coefficient. The clearest evidence of 
the enhancement of this coefficient by the presence of an Efimov trimer at zero-energy is the peak found around 
895 G by Williams et al. [3]. It is indeed located in a region of very large scattering lengths where the universal 
theory should be valid 

To calculate the recombination coefficient, we proceed as follows. We distinguish between two types of dimers: 
dimers which are included in the theory (through the solutions of l/a(in) = k), referred to as shallow dimers, 
and dimers which are not included in the theory, referred to as deep dimers. Recombination to shallow dimer 
(jk) appears in the 3-body wave function as an outgoing wave between dimer (jk) and atom i. This means that 
Xi or 0i can be written as: 



UP) = 



P 2 - Q? - is 



HP) 



P 2 -Q 2 -ie' 



(22) 



where Qi = J |(p-_E + ft?) is the relative momentum between dimer (jk) and atom i, and Mi 



J (2-rt 



P_ 

(2^)3 



= ^and^= / 



d 3 p 
(2tt)3 



Hp) 



V are factors ensuring that the dimer wavefunction 



is properly normalised to unity, or the probability V to be in the open channel, respectively. The recombination 
coefficient -fQ ec to shallow dimers (jk) is then obtained by calculating the flux of that outgoing wave [T51 12"5] : 

'^imQi)?- (23) 



K l — 
-^rec — 



Recombination to deep dimers occurs at distances on the order of the deep dimers' size, typically given by 
the range of the interactions. Therefore, the coefficient for recombination to deep dimers can be estimated by 
calculating the probabibility of finding the three atoms within that range [221 130] : 

( dee P)-£_A_/ d 3 Rd 3 r\V(R,r)\ 2 (24) 



rec 
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The results are not very sensitive to the precise choice of the range IZg, and the constant factor £ is expected to 

( deep ) 

be on the order of unity. Typically, if IZo varies by 10%, the rate Ky ec changes by 6%. As a typical size, we 

chose 7Zq = a. The total recombination coefficient is K TCC = £\ K\ cc + A'^f^. Alternatively, in the contact 
interaction model, one can set a complex three-body parameter A = |A|e I?7 , where r\ > phenomenologically 
reproduces short-distance losses due to recombination to deep dimers. In that case, the total recombination 
coefficient at zero energy is given by (see Appendix C) 

A reC = — ImVxi(O). (25) 

i 

This method also has the advantage of taking into account the broadening effect of the loss strength 77 on the 
recombination profile as a function of magnetic field. We checked that the two methods give similar results. 
The calculated recombination coefficients for different models are represented in Fig. |4| 

In the universal model, we have to adjust |A| to 1.165 a -1 and r\ to 0.016 (£ ps 0.8) in order to reproduce the 
895 G peak. This peaks correspond to the appearance of an Efimov trimer at the 3-body collisional threshold 
(zero energy) as represented in Fig. [T] Then we obtain the variation of the recombination coefficient around 895 G 
in the whole universal region B, as first calculated by Eric Braaten et al. [31]. We extend the calculation to 
low magnetic field in universal region A where earlier measurements of recombination were performed |12l 113) . 
Interestingly, the same three-body parameter roughly reproduces the variation of the recombination coefficient: 
a plateau with two peaks on both ends. However the two peaks are not exactly at the right locations. This is 
not unexpected, since different universal regions separated by zero crossings of the scattering lengths are known 
to have different 3-body phases in general [35]. Changing the value of |A| to 1.076 a -1 , and 77 to 0.115 gives a 
fair agreement between the measured recombination and the universal model - see the inset in the top panel of 
Fig. [4j Note that the two peaks correspond again to zero-energy crossings of a trimer in Fig. [l] This was in fact 
the first indication of the underlying Efimov physics in lithium 6 [33 - 35 . However the calculated recombination 
coefficient shows two marked peaks, while in the experimental data the left peak is much more pronounced than 
the right one. Variations of the loss parameter 77 with magnetic field were subsequently proposed to improve 
the agreement of the universal model with the data [36l [37] . The reason is that in the universal model the 
least-bound dimers cease to exist at low magnetic field because the scattering length becomes negative. Since 
in reality those dimers are still present, the loss parameter 77 effectively accounts for transitions to these dimers. 
As their energy varies with magnetic field, it was inferred that 77 should also change with magnetic field [37j . 

In the single-channel contact model and the two-channel separable model, the least-bound dimers are explicitly 
included over the full range of magnetic field. For binding energies larger than the typical van der Waals energy 
associated with the range the interaction, corresponding to magnetic field smaller than 600 G, their description 
becomes unrealistic. However their energy remains accurate, and they still provide a simple model for the 
recombination mechanism. In the single-channel contact model, we adjust |A| to 1.455 a -1 and 77 = 0.0033 
(£ fa 0.8) to reproduce the 895 G peak. Note that because the scattering lengths are now energy-dependent, the 
values of |A| and 77 have been altered in order to obtain the same physical situation. This is because the choice of 
the high-energy (i.e. short-distance) behaviour of the two-body interactions affects the choice of the three-body 
phase. The calculated recombination rate in the universal region B is nearly identical to the universal model. 
It is however numerically difficult to extend the calculation to the universal region A, presumably because in 
that region the dimers' binding momenta Ki are close or even exceed the cutoff momentum A. This problem 
does not occur for the two-channel separable model. 

In the two-channel separable model, we first adjusted b to 1.31 a to match the effective range of the two- 
body interaction - see Fig. [2j Suprisingly, this choice perfectly reproduces the location of the peak at 895 G. 
Again, the results in universal region B are very similar to those of the universal model. In universal region A, 
however, the results are different and significantly off the observed peak locations. Changing b to 0.75 a (thus 
using a wrong effective range, while preserving the correct scattering lengths and dimer binding energies) we 
can effectively change the 3-body phase and get better agreement - see bottom panel inset in Fig. [4] These 
numerical results also seem consistent with the semi-analytical approach of Ref. [37] based on the universal 
theory and a magnetic-field-dependent 7/. 

From all these results, we conclude that the 3-body physics of lithium 6 at zero energy is essentially consistent 
with the universal theory. 

B. Experiments at negative energy 

From the previous analysis, it is possible to predict the energy spectrum of trimers using the previously 
adjusted parameters. The spectrum based on the universal model, first predicted in Refs. [31, 34. 35 , is shown 
in Fig. [T] In universal region A, it predicts a trimer state which connects to no dimer but dissociates into three 
atoms at two magnetic field values corresponding to the two 3-body recombination peaks. In universal region 
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Figure 4: Recombination coefficient -Krec as a function of magnetic field. The stars, dots and circles indicate the 
experimental measurements of Refs. [13] - |12| . and [3], respectively. Top panel: calculations from the universal model 
for |A| = 1.165 a -1 and rj = 0.016. The inset shows results in universal region A for |A| = 1.076 a" 1 and r\ — 0.115. 
Bottom panel: calculations from the two-channel separable model for b — 1.31 a and £ = 0.08. The inset shows results in 
universal region A for b — 0.75 a and £ = 0.5. In both panels, thick black curves show the total recombination coefficient, 
dashed grey curves show the recombination coefficient to deep dimers, and coloured curves show the recombination to 
shallow dimers: (13) green, (23) blue, and (12) red. 



B, it predicts the existence of two trimers which connect to the dimer (23) at 598 G and 672 G. In fact, the 
first meeting point is outside the universal region where the dimer binding energy clearly deviates from the 
universal formula h 2 /ma 2 . Therefore, it was expected in Ref. [3T] to be unreliable. On the other hand, the 
672 G prediction is right in the universal region and was thought to be reliable. It turns out that the dimer 
energy still has a small but appreciable deviation from universality at that magnetic field, due to two-body 
finite-range corrections. However, computing the trimer energy with the single-channel contact or two-channel 
separable model (both of which include those two-body finite range corrections) leads to a similar trimer energy 
curve which again meets the dimer curve at around 672 G - see Fig. [5j In other words, the two-body finite 
range corrections shift both the dimer and trimer energies, but do not modify the magnetic field of their meeting 
point. 



1. Atom-dimer relaxation 



To check those predictions, experimentalists prepared a mixture of dimers (23) and atoms 1, and observed 
the rate of relaxation to deep dimers as atoms and dimers collide. The density n\ of atoms in state 1 therefore 
decreases according to the rate equation: 



■hi = -^ re i"-i^23, 



(26) 
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Figure 5: Trimers and dimers connection regions. Left: ground-state Efimov trimer. Right: first excited-state Efimov 
trimer. Dimer curves are the same as in Fig. n] Trimer curves are obtained for different models: universal model (dashed 
grey), single-channel contact model (purplejf two-channel separable model (pink), and single-channel model with the 
energy-dependent parameter A given in Fig. [8] (orange). All models feature a zero-energy resonance at 895 G. 



where 77-23 is the density of dimers (23) and K Ye y IS the relaxation coefficient. Calculating the relaxation 
coefficient from the theory is similar to the case of the recombination coefficient. For the case of dimer (23) 

colliding with atom 1 at energy E — Eyi — | ^ 1 1 the quantities Xi can be written as: 

Xi{P) ■ ~ 



4ttMP) 
P 2 - Q? - is 



M 6 tl (2ir) i 5 i (P-Q 1 ) + 



(27) 



corresponding to an incident wave and outgoing waves of amplitude /j (Qi) and momentum Qi — W | ( p- E + re?) . 

The relaxation coefficient to shallow dimers (jk) is then obtained by calculating the flux of the corre- 

sponding outgoing wave: 



K 



Ah 



rel 



i\fi(Qi)f 



and the total relaxation coefficient is 



4/i 



K Yel = — (lm/i(Qi) " Qi\h(Qi)\ 2 



(28) 



(29) 



Similarly to 3-body recombination, whenever the magnetic field is close to the meeting point between a dimer 
and trimer, a resonance occurs which strongly enhances the relaxation coefficient. Indeed, two peaks were 
observed near the expected meeting points, at 602 G and 685 G [7J - see Fig. [6] However, the significant 
deviation of the second peak location from the expected value 672 G is somewhat surprising. Indeed, we 
checked that if we change the 3-body parameter to obtain a peak at the measured 685 G, then the peak in 
the 3-recombination coefficient moves to 865 G, which seems incompatible with the measured value of 895 G in 
Ref. 0. 

Similar measurements of the relaxation coefficient involving a dimer (12) colliding with atom 3 [7] revealed 
some dips in the relaxation coefficient at 610 G and 690 G - see Fig. [7| These dips are expected to result 
from two-pathway interferences related to the Efimov physics of the system |llj . Actually only the first dip 
appears as a clear signature of such interference. The second dip, as our calculations suggest, may be due to the 
combined effect of relaxation to shallow dimers and deep dimers. While both the universal and single-channel 
contact models predict a dip at 600 G, the real dip is again shifted away from the theoretical expectation. 

Since all non-universal two-body corrections have been taken into account in the theoretical models, we 
concluded in Rcf. [5] that a non-universal three-body correction is also needed. In the contact model, non- 
universal two-body corrections arise from the energy dependence of the scattering lengths (p) . In the same 
fashion, we expect the three-body parameter A to be energy-dependent in general. By adjusting the value of 
A(p) to fit each peak and dip locations we obtain some insight on this energy dependence. The result is plotted 
in Fig. [8] and looks consistent with a smooth but non-linear variation with energy. By interpolating A(p) as a 
function of p, we can recalculate all previous curves and obtain reasonable agreement with experimental data, 
see Figs.[TJ[6j and[7j This adjusted non-universal 3-body model predicts that the ground-state Efimov trimer 
at large magnetic field is shifted by about 20 MHz from the universal prediction. Direct measurement of that 
energy would clearly validate or invalidate our assumption of the energy dependence of the 3-body parameter 
A. The fact that the 3-body parameter of an effective theory is energy-dependent to describe the ground-state 
trimer is not suprising, since ground states of Efimov series are always non-universal. It is less obvious, however, 
why it is also energy-dependent near the second trimer, which should be closer to universality. 
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Figure 6: Relaxation coefficient for dimer 23 colliding with atom 1 as a function of magnetic field. The grey and blue 
dots indicate the measurements from Refs. [7] and [6]. Left: universal model with |A| — 1.165 a -1 and r/ = 0.016 - 
see similar calculations in Refs. [3 1 1 138| . Middle: contact model with |A| = 1.455 a -1 and r] = 0.0033. Right: contact 
model with the energy-dependent |A| of Fig. [8] and r\ — 0.0033 (dashed curves: rj — 0.0400). Black curves show the total 
relaxation coefficient, red curves show the relaxation to dimer 12, and grey curves show the relaxation to deep dimers. 
The horizontal line indicates the unitarity limit given by the typical collisional energy (~100 nK) in the experiments. 




575 600 625 650 675 700 725 750 575 600 625 650 675 700 725 750 575 600 625 650 675 700 725 750 
Magnectic field (gauss) Magnectic field (gauss) Magnectic field (gauss) 



Figure 7: Relaxation coefficient for dimer 12 colliding with atom 3 as a function of magnetic field. The red dots show 
the results of the measurements from Ref. [7] . The same conventions as in Fig. [6] are used. 



2. Atom-dimer association spectroscopy 



Very recently, the binding energy of the excited Efimov trimer was directly observed by association spec- 
troscopy |15[l28j. The experimentalists prepared a mixture of dimers and atoms, and applied a radio-frequency 
field to induce a transition to the trimer state. The preliminary results [15] reported excellent agreement with 
our own theoretical prediction [6 of the binding energy based on the single-channel contact model with the 
energy-dependent three-body parameter A(p) - see Fig. [9] However, it turns out that the thermal shift of the 
association peaks is non- negligible in this experiment, and was not taken into account in the reported measure- 
ments of the binding energies. Correcting for these shifts makes the measurements deviate significantly from the 
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Figure 8: Energy dependence of the three-body parameter A (in units of a^ 1 ) of the single-channel contact model adjusted 
to fit measured resonances. Each dot corresponds to an adjustment to a peak or dip of the relaxation coefficient. 
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Figure 9: Direct measurements of the second trimer energy: dots correspond to measurements of Ref. [15) . which were 
taken at about 1 \iK , and squares show the results of the measurements in Ref. 28 which were taken at lower temperature 
where thermal shifts are negligible. Using the same conventions as Fig. [I] thick curves represent the dimer energies and 
the orange curve represents the trimer energy obtained from the single-channel contact model with the energy-dependent 
3-body parameter shown in Fig. [8] 

original prediction. Worse still, to account for the deviation, the parameter A(p) should be adjusted in a way 
which breaks its expected smooth variation with energy. This problem is the subject of a separate study [28] . 

IV. CONCLUSION 

In this paper we provided an overview of the various experimental results concerning the Efimov physics in 
a three-component lithium 6 system. Using different theoretical models, we found that the Efimov features 
measured at nearly zero energy are essentially consistent with the so-called Efimov scenario based on the uni- 
versal theory developed in Vitaly Efimov's original papers. However, the features measured at negative energy, 
although qualitatively consistent with the Efimov scenario, show some significant deviations from universal 
theory, as well as theories which fully take into account non- universal two-body corrections. To account for 
this, we introduced phenomenological non-universal 3-body corrections through an energy dependence of the 
three-body parameter adjusted the experimental data. While most measurements can be explained by this ad 
hoc parametrization, it does not seem to be consistent with the very recent measurements of a trimer's binding 
energy, and further studies are needed to fully undertand the physics at negative energy. The two-body model 
which the present study is based on is thought to be very accurate. However, even minor inaccuracies at the 
two-body level has been shown to be important in the three-body physics of other systems [39], and therefore 
deserve further investigation in our system as well. 



Appendix A 



Here, we give the derivation of Eq. (10). From Eq. ([9]), we perform the inverse Fourier transformation with 
respect to the second variable: 



4f(P,r f ) = 4f (P,i f ) 



-ip\A 



-Xi(P) + 4tt 



P 2 + p 2 - mE/h 2 + is (2vr) 3 



(30) 



Then, applying the operator lim r _j.o 75>( r ')j an d using the definition of xi, we find: 



ai(ryp) 



UP) 



r->0 



X 2 (-P 2 )+X3(J 3 3 ) d 3 p 
\P 2 +p 2 - mE/h 2 + ie (2tt) 3 



(31) 



Using P2 = — 7t — p, and P3 = — + p, we can change the integration variable to get 



Xi(P) = ^ (P,r 



r-f0 



-1pXi(P) +4tt 



X2(q) + X3(q) 



d 3 q 



IP 2 + (q+ f ) 2 - mE/h 2 + ie ( 2 ^) 3 



(32) 
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IP 



Xi(P)-4tt 



X2(q) + X3{q) 



d 3 q 



P 2 +q 2 + P-q- mE/h 2 + ie (2tt) 



r->0 



(33) 



Assuming that Xi(9) — Xi(<?)j an d writing the integration element d 3 q = q 2 dqdipd{— cos 9), with 
P ■ q = Pq cos 9, we can perform the integration over the angles (p and 0: 



1 



7p 



47T 

8^ 



ai(«7p) 
d^p 



Xi(P) 



d{-cos9) 



1 P 2 + q 2 - Pq cos 6» - mE/h 2 + ie 



(X2(q)+xM)Q 2 dq = *o(P,^ 



r-s-O 



(34) 



^—,+ip)xi(p)-- [ ^^ P l + S + P p mE JZ fe(g) + x 3 (g))g 2 ^ = \MP, rl] , (35) 

«i( ? 7p) / n J Pq P 2 + q z — Pq — mE /h 2 L J r-yo 



from which we obtain Eq. ( 10 ) after introducing a cutoff A to the integral. 



Appendix B 



Here we derive the connection between the upper bound A of the integral in Eq. ( 10 ) and the boundary 
condition at short distance between the three particles. According to Danilov [20] and Min os-Faddeev [21], the 
solutions Xi(p) of the original Skorniakov - Ter-Martirosian equations (without upper bound A) are superposi- 
tions of two linearly independent solutions, which have the following asymptotic form: 



XiiP) > A 

p—too 



sin(s lnp) cos(s lnp) sin(s ln^) 



B 



oc 



P 



P 



P 



(36) 



where sq m 1.00624 is a constant. Thus, some extra condition can determine one particular solution (up 
to normalisation) by fixing the ratio B/A — tan(so hi Ao). The quantity Ao is the 3-body parameter of the 
original Efimov theory. It is set by imposing a boundary condition on the wave function at short distance pQ. 
Alternatively, one can cut off the integral in Eq. ( 10 ) at some high momentum A [55J [5TJ [35] , and this imposes 

be zero: 



the condition that the remaining part of the intcgra 



10 q P 2 + q 2 + Pq- mE/h 2 
q P n P 2 + q 2 -Pq- mE/h 2 



(x 3 (q) + xk{q)) 



sin(s ln^ 

k 2 



-dk = 



This condition selects a particular solution. By evaluating the above integral, one obtains 



s cos(s In + sin(s In 







(A/Ao)(l + *g) 

which gives an explicit relation between the three-body parameter Ao and the imposed cutoff A: 

arctan sq + im ... 

Ao = Aexp( ), with n integer. 

so 



(37) 



(38) 



(39) 



Appendix C 



The expression for the total loss rate coefficient at zero energy, Eq. (25), can be seen as a consequence of 
the optical theorem. Here we indicate a short derivation. We start from the wave function in Eq. |9]) in space 
coordinates at zero energy 



d 3 P 3 d 3 pj 



4tt 



(2tt)3 (2vr)3 \P 2 + p 2 + ie 



UPj) 



(40) 



For large Rj, i.e. when all three atoms are far from each other, it can be approximated as ^(i?, r) = 1 
GoiR^EjXM, where 



G Q (R,f) = 



d 3 Pj d 3 pj 



4n 



(2tt) 3 (2ir) 3 ipj+fl + ie 



e iP 3 -Rj e ipj-rj _ 



2y/3 1 

(3ir/2) 2 W 



(41) 
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and TZ — yr 2 + |i? 2 is the hyper-radius of the system. Calculating the flux of probability current of the wave 
function through a large hypersphere S, we obtain 



2h Im f vTV TC * • dS = — 



3-o5 



2H ( ( V3 Y dG ° 2n n 

~ \ "a" ) ~diz r[3] 



Ah 



= -^>5>(0). 

Til < J 



Im^ Xl (0) 



(42) 
(43) 
(44) 



This incoming flux should balance the outgoing fluxes in other sectors where two atoms are recombined into 
a shallow dimer [Eq. (23)], as well as the loss induced by the imaginary part of the three-body parameter A, 
which phenomenologically describes recombination to deep dimers. Therefore, Eq. (44) represents the total 
recombination coefficient K^ec- 
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